o 

(N 



> 

On 



o 



Recent progress in the theory of homogenization with 
oscillating Dirichlet data 

David Gerard- Varet*and Nader Masmoudi^ 
January 31, 2013 



In this talk we study the homogenization of elliptic systems with Dirichlet boundary 
condition, when both the coefficients and the boundary datum are oscillating, namely e- 
periodic. In particular, in the paper [9], we showed that, as e — > 0, the solutions converge in 1? 
with a power rate in e, and we identified the homogenized limit system and the homogenized 
boundary data. Due to a boundary layer phenomenon, this homogenized system depends in 
a non trivial way on the boundary. The analysis in [U] answers a longstanding open problem, 

■ raised for instance in [3]. 
<' 

■ 1 Introduction 

Homogenization of elliptic systems arises in several physical problems where a mixture is 
present. Some of the main applications of the theory are the diffusion of heat or electricity 
in a non- homogeneous media, the theory of elasticity of mixtures, ... Physically, the main 
goal of the theory is to try to compute accurate and effective properties of these mixtures. 
Mathematically, we have to find a limit system towards which the solutions of homogenization 
problem converge. This passage from "microscopic" to "macroscopic" description is called in 
the literature "homogenization". 

When both the coefficients of the system and the boundary datum are oscillating (e- 
(T) ■ periodic) and due to a boundary layer phenomenon, this homogenized system depends in a 

non trivial way on the boundary. In this talk, we answer a longstanding open problem, raised 
for instance by Bensoussan, Lions and Papanicolaou in their book "Asymptotic analysis for 
^ ■ periodic structures" [4, page xiii]: 

5—( < 

Of particular importance is the analysis of the behavior of solutions near bound- 
aries and, possibly, any associated boundary layers. Relatively little seems to be 
known about this problem. 

In particular this result extends substantially previous works obtained for polygonal do- 
mains with sides of rational slopes as well as our previous paper [8] where the case of irrational 
slopes was considered. We hope that these notes give a better understanding of the proof of 
the result in |9|. 
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2 The homogenization problem 



We consider the homogenization of elliptic systems in divergence form 



V • (A (-/e) Vit) (x) = f, 



(2.1) 



set in a bounded domain Q, of M. d , d>2, with an oscillating Dirichlet data 



u(x) = ip(x, xje) 



x G on. 



(2.2) 



As is customary, e > is a small parameter, and A = A{y) takes values in M^(Mjv(M)), 
namely A afS (y) G M N (R) is a family of functions of y G K d , indexed by 1 < a, (3 < d, with 
values in the set of N x N matrices. Here, u = u(x) and tp = <p(x, y) take their values in R . 
We recall, using Einstein convention for summation, that for each 1 < i < N, 



In the sequel, Greek letters a, f3, ... will range between 1 and d and Latin letters i, j, k, ... will 
range between 1 and N. 

In the context of thermics, d = 2 or 3, iV = 1, u is the temperature, and a = A(-/e)\7u is 
the heat flux given by Fourier law. The parameter e models heterogeneity, that is short-length 
variations of the material conducting properties. The boundary term ip in (|2.2p corresponds 
to a prescribed temperature at the surface of the body and / is a source term. In the context 
of linear elasticity, d = 2 or 3, N = d, u is the unknown displacement, / is the external load 
and A is a fourth order tensor that models Hooke's law. 

We make three hypotheses: 

i) Ellipticity: For some A > 0, for all family of vectors f = £f G R Nd 



ii) Periodicity: Vy G Mr, V/i G Z , Vx G dU, A(y + h) = A(y), ip{x,y) = ip(x,y + h). 

iii) Smoothness: The functions A, f and 99, as well as the domain £1 are smooth. It is actually 

enough to assume that (p and are in some H s for s big enough, but we will not try 
to compute the optimal regularity. 

The main question we are trying to answer is the following: 

Question: What is the limit behavior of the solutions u £ as e — > ? Can we go beyond 
the limit and compute a full expansion of u e ? 

This question goes back at least to the 1970's, and a classical approach consists in trying 
a two-scale expansion: 

Classical approach: Two-scale asymptotic expansion: 



{V-A(./e)VuUx) :=d Xa A«f (./e) d x „ Uj (x). 




a 



a 



app 



u°(x) + eu (x,x/e) + ... + e n u n (x,x/e) 



(2.3) 



with u l = u l {x,y) periodic in y. 
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3 Case without boundary 



The two-scale approach works well in the case without boundary, namely in the whole space 
case or in the case of a periodic domain (say of period 1 and e is taken to be equal to 1/n with 
n an integer). In particular one can construct inductively all the terms in the expansions. 
Let us recall few classical facts (see for instance [231 EH H3 E]) : 

i) The construction of the u l, s involves the famous cell problem 



-V • {AVx 1 ) (y) = V a ■ A ai (y), y in T d (3.1 



with solution x 1 6 M N (R). 

ii) The solvability condition for u 2 yields the equation satisfied by u , namely u° (which does 
not depend on y) satisfies 

V-A°Vu° = f (3.2) 
where the constant homogenized matrix is given by 



A o, a p = / A *fi( y)dy+ / A^(y)d y y(y)dy. 



The second term in the expansion (|2.3|) reads 

u\x,y) := u 1 (x,y)+u 1 (x) := - X a (y)d Xa u°(x) + u\x), (3.3) 



where x is again the solution of (|3,ip . 

To find an equation for the average part n 1 (x), one needs to introduce another family of 
1-periodic matrices 

^ = ^( s )eI n (I),^ = l,.,(i, 

satisfying 

- V,, • AV y T a P = B a P - I B a ^ , f T a/3 = 0, (3.4) 

where 

£>a/3 ._ _ ^07 _ ^ 

Formal considerations yield 

u 2 (x,y) := T a ^^- - X a d a u l +u 2 (3.5) 

and that the average term u 1 = v}(x) formally satisfies the equation 

- V • A^u 1 = c a ^ - - — , c a ^ := / A 171 -— A^x 1 - (3.6) 

dx a dxi 3 dx 1 J y dy v 

We refer to [2] for more details. 

Inductively, one can keep constructing all the terms of the expansion by introducing new 
corrector families as in (|3.4p and solving homogenized systems to determine u k as in 
Note that in this case, we do not need an extra boundary condition to solve (|3.6p . 
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4 Case with boundary 

Two boundary conditions have been widely studied and are by now well understood as long 
as we are only interested in the first term of the expansion: 

1. The non-oscillating Dirichlet problem, that is (|2.ip and (|2.2p with <p = (p{x). 

2. The oscillating Neumann problem, that is (|2.ip and 

n(x) ■ (A(-/e)Vu) (x) = (p(x,x/e), x G dfl, (4.1) 

where n(x) is the normal vector and with a standard compatibility condition on ip. Note 
that in thermics, this boundary condition corresponds to a given heat flux at the solid 
surface. 

Notice that in both problems, the usual energy estimate provides a uniform bound on the 
solution u £ in H l {Vi). 

For the non-oscillating Dirichlet problem, one shows that u £ weakly converges in H l (£l) 
to the solution u° of the homogenized system 

f -V • (A°Vu°) (x) = f, xGQ, 

\ n ( 4 - 2 ) 

[ u (x) = (p(x), x G d£l. 

It is also proved in [4] that 

u z( x ) = u °(x) + su 1 {x,xle) + 0(y/e), in H\Q). (4.3) 

Actually, an open problem in this area was to compute the next term in the expansion in the 
presence of a boundary, namely to compute n 1 (x,x/e). Indeed, it is not difficult to see that 

u\x,y) = - x a (y)d Xa u°(x)+u 1 (x), (4.4) 

where v}(x) solves the homogenized equation (|3.6p . However, the main difficulty is to find 
the boundary data for ■u 1 (x). The new analysis of [S] gives an answer to this problem (see 
also next section). 

For the oscillating Neumann problem, two cases must be distinguished. On one hand, if 
(90 does not contain flat pieces, or if it contains finitely many flat pieces whose normal vectors 
do not belong to WL n , then 



ip(-, -/e) — > ip := / (p weakly in L 2 (<9$7) 
J\o,i] d 



'[0,1] 

and u e converges weakly to the solution u° of 

f -V • (A°Vu°) fx) = 0, 

\ nix) ■ (^°Vn°) (x) = Tp(x), x G dVL. 

On the other hand, if dQ does contain a flat piece whose normal vector belongs to IRQ^, then 
the family </?(•, -/e) may have a continuum of accumulation points as e — > 0. Hence, u £ may 
have a continuum of accumulation points in H 1 weak, corresponding to different Neumann 
boundary data. We refer to [4] for all details. 
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5 Case of an oscillating Dirichlet data 

Here we study (|2.1|) with the boundary data (|2.2|) . One of the motivation to study this case 
is actually to understand the boundary condition for u x (x) which appears in (|3.3j) . 

Let us explain the two main sources of difficulties in studying ()2.1() - (|2.2|) : 

i) One has uniform LP bounds on the solutions u £ of (|2.ip - (|2.2p . but no uniform H 1 bound 

a priori. This is due to the fact that 

\\x (->• <p(x,x/e) Htfi/^an) = 0(e _1/2 ), resp. \\x i-> v(x,x / e)\\ LP{m) = 0(1), p > 1. 
The usual energy inequality, resp. the estimates in article |3J page 8, Thm 3] yields 
\\u £ \\m{n) = 0(e~ 1/2 ), resp. \\u £ \\ LP{n) = 0(1), p > 1. 

This indicates that singularities of u £ are a priori stronger than in the usual situations. 
It is rigorously established in the core of the paper [9]. 

ii) Furthermore, one can not expect these stronger singularities to be periodic oscillations. 

Indeed, the oscillations of cp are at the boundary, along which they do not have any pe- 
riodicity property. Hence, it is reasonable that u £ should exhibit concentration near dtl, 
with no periodic character, as e — > 0. This is a so-called boundary layer phenomenon. 
The key point is to describe this boundary layer, and its effect on the possible weak 
limits of u £ . 

It is important to note that there is also a boundary layer in the non-oscillating Dirichlet 
problem, although it has in this case a lower amplitude (it is only necessary to compute the 
boundary data of u 1 to solve (|3.6p ). More precisely, it is responsible for the 0{yfe) loss in 
the error estimate f|4.3j) . If either the L? norm, or the H 1 norm in a relatively compact 
subset u <s 0, is considered, one may avoid this loss as strong gradients near the boundary are 
filtered out. Following Allaire and Amar (see [2J Theorem 2.3]), we can give a more precise 
description than (|4.3p : 

u £ = u°( x ) + 0(e) in L 2 (n), u £ (x) = u°(x) + eu 1 ^ x/e) + 0(e) in H 1 (oj). (5.1) 

Still following [2], another way to put the emphasis on the boundary layer is to introduce the 
solution ulf(x) of 

' -V • A (-) Vulf = 0, ieS]cR d , 

W bl (5.2) 

ulf = — n 1 (x, x/e), x G 00. 

Actually, understanding this system and requiring that u^f goes to zero inside the domain f2 
allows to determine the right boundary condition for u . Hence, one can show that 

u £ (x) = u °(x) + eu 1 (x,x/e)+eu 1 b f(x) + O(e), in H 1 ^). (5.3) 

or 

u £ (x) = u°(x)+eu 1 (x,x/e)+eu 1 b f(x) + 0(e 2 ), in L 2 (fl). (5.4) 

Note that system (|5.2p is a special case of (|2.1|) - (|2.2p . Thus, the homogenization of the 
oscillating Dirichlet problem may give a refined description of the non-oscillating one. 



5 



6 Prior results 



Until recently, results were all limited to convex polygons with rational normals. This means 
that 

U := nfLi l x , n k ■ x > c fc | 

is bounded by K hyperplanes, whose unit normal vectors n k belong to RQ d . Under this 
assumption, the study of (|2.ip - (|2.2p can be carried out. The keypoint is the addition of 
boundary layer correctors to the formal two-scale expansion: 



u £ (x) ~ u°(x) + £U 1 (x,x/e) + y~]vg (x,-J , (6.1) 

k £ 

where vt, = v k ,(x,y) G M n is defined for and y in the half-space 



n 



e,k 



{y, n k ■ y > c k /e} 



These correctors satisfy 



(-Vy A(y)V y v k bl = 0, ye^, 
\ vk = <p(x,y)-u°(x), y£dtf> k . 

We refer to the papers by Moskow and Vogelius [19], and Allaire and Amar [2 J for more 
details. These papers deal with the special case (|5.2p . but the results adapt to more general 
oscillating data. Note that x is just a parameter in (|6.2p and that the assumption n k G 
M.7j d yields periodicity of the function A(y) tangentially to the hyperplanes. The periodicity 
property is used in a crucial way in the aforementioned references. First, it yields easily 
well-posedness of the boundary layer systems (|6.2p . Second, as was shown by Tartar in [18, 
Lemma 10.1] (see also subsection I7.2p . the solution t^(x,y) converges exponentially fast to 
some vt, ^(x) = <p k (x) — u°(x), when y goes to infinity transversely to the k-th hyperplane. 
In order for the boundary layer correctors to vanish at infinity (and to be o(l) in L 2 ), one 
must have vh ^ = 0, which provides the boundary condition for vP. Hence, u° should satisfy 
a system of the type 

(-V-(A°Vu°)(x) = f, xefi, 
I u°(x) = <p#(x), x G 8ft. 

where <p*(x) := tp k (x) on the k-th side of Q. Nevertheless, this picture is not completely 
correct. Indeed, there is still a priori a dependence of (p k on e, through the domain Q e,k . In 
fact, Moskow and Vogelius exhibit examples for which there is an infinity of accumulation 
points for the ip k, s, as e — > 0. Eventually, they show that the accumulation points of u e in L 2 
are the solutions u° of systems like (|6.3p , in which the ip k 's are replaced by their accumulation 
points. See [19J for rigorous statements and proofs. We stress that their analysis relies heavily 
on the special shape of Q, especially the rationality assumption. 

A step towards more generality has been made in our recent paper [8] (see also [7]), in 
which generic convex polygonal domains are considered. Indeed, we assume in [8] that the 
normals n = n k satisfy the Diophantine condition: 

For all ^Z d \ {0} \P n ±(£)\ > for some k, I > 0, (6.4) 



6 



where P n ± is the projector orthogonal to n. Note that for dimension d = 2 this condition 
amounts to: 

For all ££Z d \ {0} {n 1 - ■ f | := | - n 2 £i + ni&| > « If |"', for some «, Z > 0, 
whereas for d = 3, it is equivalent to: 

For all f G Z d \ {0} |n x f | > « |f |~', for some k, I > 0. 

Condition (|6.4p is generic in the sense that it holds for almost every n G S 1 . 

Under this Diophantine assumption, one can perform the homogenization of problem (|2.1|) - 
()2.2p . Stricto sensu, only the case (j5.2jl . cZ = 2,3 is treated in [8], but our analysis extends 
straightforwardly to the general setting. Despite a loss of periodicity in the tangential variable, 
we manage to solve the boundary layer equations, and prove convergence of away from the 
boundary. The main idea is to work with quasi-periodic functions instead of periodic ones 
(see also subsection 17. 3p . Interestingly, and contrary to the "rational case", the field <p* does 
not depend on e. As a result, we establish convergence of the whole sequence u e to the single 
solution vP of (|6.3|) . We stress that, even in this polygonal setting, the boundary datum ip* 
depends in a non trivial way on the boundary. In particular, it is not simply the average of 
ip with respect to y, contrary to what happens in the Neumann case. 

7 Main new result and sketch of proof 

The main new result of [9] is to treat the case of a smooth domain: 
Theorem 1 (Homogenization in smooth domains) 

Let Q be a smooth bounded domain o/R d , d>2. We assume that it is uniformly convex (all 
the principal curvatures are bounded from below). 

Let u £ be the solution of system (|2.ip - (|2.2p . under the ellipticity, periodicity and smoothness 
conditions i)-iii). 

There exists a boundary term tp* (depending on ip, A and £1), with ip* G L p (dQ) for all finite 
p, and a solution u° of (|6.3p . with u° G L P (Q) for all finite p, such that: 

IK - < C Q e a , for all 0<a< ^=1. (7.1) 

We will present a sketch of the proof of theorem [1] 

^From the two difficulties explain in section [5l we know that the first term in the expansion 
(I2.3P should be independent of y and should solve (|3.2p . The main question is : 
Question: What is the boundary value <p° of u° ? 
Solution : We need a boundary layer corrector 
Difficulty: There is no clear structure for the boundary layer. 
Guess : The boundary layer has typical scale e and there are no curvature effect: 

• Near a point xq G d£l, we replace d£l by the tangent plane at xq: 

To(dQ) := {x, x ■ uq = xq ■ no} 
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We dilate by a factor e . 



Formally, for x ~ xq, one looks for 



u £M {x) « C/ (x/e) 



where the profile Uq = Uo(y) is defined in the half plane 



{y, y-n > e 1 x • n }. 



It satisfies the system: 



V„ ■ (AVyU ) = in fig, 

Uo\dH$ =l f- <P°{xo)- 



Notice that in this system, xo is just a parameter. 

7.1 Study of an auxiliary boundary layer system 

The previous heuristic justifies the study of 




where H := {y, y ■ n > a} and (j) is 1-periodic in y. 
We expect that the solution U of ()BL|) satisfies: 



U^U n 



as y • n — )• +00, 



(7.2) 



(BL) 



for some constant Uoo = Uoo(<p) that depends linearly on (j). 

If we go back to Uq which solves (|7.2p . one can derive the homogenized boundary data 
<^°. Indeed: 

• On one hand, one wants Uq — >• (localization property) when y ■ n +00. 

• On the other hand, 

->■ C/oo(</? - <p°(x )) = Uoo((p) - (p°(x ) 

so that we need to take: 

^ (x ) := Uoo{y). 
This formal reasoning raises many problems : 
1. The well-posedness of (jBLp is unclear: 

- No natural functional setting (no decay along the boundary). 

- No Poincare inequality. 

- No maximum principle. 
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2. The existence of a limit Uoq for (|BL|) is unclear: 
There is an underlying problem of ergodicity. 

3. Uqo depends also on H, that is on n and a: 

- There is no obvious regularity of 11^ with respect to n. 

- Back to the original problem, our definition of cp (xo) depends on xq, but also on the 
subsequence e. Indeed, there is possibly many accumulation points as e — > (see [TO]). 

7.2 Polygons with sides of rational slopes 

In this cases, the boundary layer systems of type (|BLp can be fully understood (see [191 E])- 
For simplicity, we only concentrate on the case d = 2. 

1. Well-posedness: The coefficients of the systems are periodic tangentially to the boundary. 
After rotation, they turn into systems of the type 



V z • {BV Z V) = 0, z 2 >a, 
V\ Z2=a = ip, 



with coefficients and boundary data that are periodic in z\ which yields a natural 
variational formulation. 

2. Existence of the limit : Saint- Venant estimates on (jBLip . 

One shows that F(t) := J Z2>t |V 2 V| 2 dz satisfies the differential inequality. 

F(t) < -CF'(t). 

/.From there, one gets exponential decay of all derivatives, and the fact that: 

V — > Vqo, exponentially fast, as z 2 — > +oo 
and hence going back to (BL), we get 

U Uoo, exponentially fast, as y ■ n — > +oo. 

A Key ingredient in this case is the Poincare inequality for functions periodic in z\ with 
zero mean. 

3. In polygonal domains, the regularity of C/qo with respect to n does not matter. However, 
for rational slopes, the limit U°° does depend on a. This means that if we go back to our 
original problem (in polygons with rational slopes), The analogue of our theorem is only 
available up to subsequences in e. Moreover, the boundary data of the homogenized 
system may depend on the subsequence. Indeed, there are examples with a continuum 
of accumulation points (see [IS]). 



(BL1) 
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7.3 More general treatment of (IBLj) 



It is worth pointing out that one can not be fully general: The existence of Uoq requires some 
ergodicity property. A simple example is : 

AU = in {y 2 > 0}, U\ V2=0 = <j> . 



If 4> l-periodic, then U(0, 2/2) — > J 4> exponentially fast. 
But there exists <j) G L°° such that U(0, 2/2) no limit. 



Indeed, we have an explicit formula: U(0, 2/2) 



1 



</2 



-</)(i) (it. For 4> with values 



7T y\ + * 2 

in {+1,-1}, the asymptotics relates to coin tossing. Hence, we need some extra structure 
(or ergodicity) to solve the problem. 

In our case, we have some ergodicity property ! For general half planes, the coefficients 
of (|BLp or (|BL1|) are not periodic, but they are quasiperiodic in the tangential variable. We 
recall that a function F = F(z\) is quasiperiodic if it reads 



F( Zl ) = T(Xzi) 



where A G and T = J- (9) is periodic over R D ( D > 1). As an example: For (jBLTjl . 
D = 2 , and A = n 1 - (the tangent vector). 

Notice that the previous results (subsection 17. 2\ correspond to the case: n £ RQ 2 . Now, 
we replace this by the small divisor assumption: 



(H) 3k > 0, |n-£| > ~\ G Z \ {0} 



Note that the assumption (H) is generic in the normal n: It is satisfied by a set of full 
measure in S 1 . But it does not include the previous result of subsection 17.21 

Proposition 2 If n satisfies (H), the system (|BLj) is "well-posed", with a smooth solution U 
that converges fast to some constant Uoo. Moreover, Uqo does not depend on a. 

Proof of the proposition: 

1. Well-posedness: involves quasiperiodicity. One has: 




where B(z) = B{\zi,Z2), "ip(z) = T(\z\,Z2). 

Notice that the functions B = B(9,t) and V = V(6,t) are periodic in 9 G T 2 . 
The idea is to consider an enlarged system in 9, t, of unknown V = V(9, t): 




(BL2) 
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where D is the "degenerate gradient" given by D = (A • Vg,dt) 
Advantage: Back to a periodic setting (6 E T 2 ). 

Drawback : We have a degenerate elliptic equation. However, we are still able to prove 
the following : 

- Variational formulation with a unique weak solution V. 

- One can prove through energy estimates than V is smooth. 

- Allows to recover V through the formula V(z) = V(Xzi, Z2). 

2. To prove the convergence to a constant at infinity, we rely again on Saint- Venant type 

estimates, adapted to (fBL2]) . Thanks to (H), we prove that F(t) := / \DV\ 2 dO dt' 

Jt>>t 

satisfies 

F(t) < C(-F'(t)) a , Va < 1. 

But, we have only polynomial convergence towards a constant. 

We point out that this better understanding of the auxiliary boundary layer systems allows 
to handle the generic polygonal domains in the next subsection. 

7.4 Extension to smooth domains 

The are at least three main difficulties to extend the previous analysis to smooth domain : 

1. The none smoothness of Uoo with respect to n. Indeed, Uoo is only defined almost 
everywhere (diophantine assumption). 

Idea : For any k > 0, we can prove that Uoo is Lipschitz when it is restricted to 

A K : = jneS 1 , |n-e| > V£ G I? \ {0}| . 

Moreover, we have that \A C K \ = 0{k). 

In the course of the proof, the construction of the boundary layer corrector can be per- 
formed in the vicinity of points x such that n(x) G A K . In some sense, the contribution 
of the remaining part of the boundary is negligible when n <C 1. More precisely, 

2. We have to approximate the smooth domains by some polygons with sides having normal 
vectors in the set A K . In doing so, we will introduce another small parameter e a . 

3. We have to construct a more accurate approximation due to the many errors made in 
the previous two points. 

Broadly, optimizing in k, a and e yields a rate of convergence. We refer to [S] for the 
details. 
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8 Conclusions 



We would like to conclude by mentioning a few related results. Recently there was many 
activity in the theory of homogenization and many new problems were addressed. We would 
like to mention some of them since we think they may give a better understand of our result 
or/and may be combined with our result: 

• Our results on the boundary data problem were recently extended to the eigenvalue 
problem, see [21J. Also, the behavior of the reduced boundary layer system (|BLj) was 
recently investigated by C. Prange in [22], without any diophantine assumption. 

• The Avellaneda-Lin type estimates were extended to the case of Neumann boundary 
conditions by Kenig, Lin and Shen [151 [151 HE] ( see a ^ so [5] for a related work). These 
estimates should be helpful to study the next order approximation for the Neumann 
boundary condition case 

• Many new probabilistic results were proved when an interface is present (see [12]) or in 
the trying to compute the accurate value of the homogenized matrix (see |1Q[ [TT] ) . 

• Some different method was used to compute homogenized boundary data for none os- 
cillating coefficient ( [TTJ [1] ) . 
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